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^ ■ Abstract 

The prediction of standard chemical evolution models of higher abundances of 3 He at 
the solar and present-day epochs than are observed indicates a possible problem with the 
yield of 3 He for stars in the range of 1-3 M . Because 3 He is one of the nuclei produced in 
J> ■ Big Bang Nucleosynthesis (BBN) , it is noted that galactic and stellar evolution uncertainties 
^ ■ necessarily relax constraints based on 3 He. We incorporate into chemical evolution models 
which include outflow, the new yields for 3 He of Boothroyd & Malaney (1995) which predict 
that low mass stars are net destroyers of 3 He. Since these yields do not account for the high 
3 He/H ratio observed in some planetary nebulae, we also consider the possibility that some 
fraction of stars in the 1-3 M Q range do not destroy their 3 He in theirpost main-sequence 
phase. We also consider the possibility that the gas expelled by stars in these mass ranges 
does not mix with the ISM instantaneously thus delaying the 3 He produced in these stars, 
according to standard yields, from reaching the ISM. In general, we find that the Galactic 
D and 3 He abundances can be fit regardless of whether the primordial D/H value is high 
(2 x 1(T 4 ) or low (2.5 x 10" 5 ). 



1 Introduction 



The abundances of the light elements, D, 3 He, 4 He, and 7 Li, serve as a critical test to the 
standard hot big bang model (Walker et al. 1991). In the case of two of these elements, 4 He 
and 7 Li, primordial values can be reasonably well determined directly from observations. 4 He 
may be inferred from low metallicity Hn regions (see e.g. Pagel et al. 1992; Olive & Steigman 
1995; Olive, Skillman, & Steigman 1996; Skillman et al. 1996; Olive & Scully 1996), while 
the uniform abundance of 7 Li in halo dwarfs is interpreted to be the primordial value for 
7 Li/H (Spite & Spite, 1982). Although possible stellar depletion of as much as a factor of 2 
can not be categorically excluded, higher depletion is in conflict with 6 Li observed in these 
stars and with maintaining the tightness of the plateau (Steigman et al. 1993; Vauclair & 
Charbonnel 1995; Lemoine et al. 1996). Fields & Olive (1996) and Fields et al. (1996) have 
argued that on the basis of these two isotopes alone, one can constrain the single parameter 
(the baryon-to-photon ratio, rj) of standard big bang nucleosynthesis (SBBN) with the degree 
of constraint depending on the allowed 7 Li depletion and the maximum systematic errors 
allowed in the 4 He determination. 

The consistency of the two remaining elements, D and 3 He, with such a best fit rj based 
on 4 He and 7 Li, may pose a challenge for many chemical evolution models. Let us first 
discuss deuterium which may eventually give the cleanest constraints but at present has 
some ambiguities. Deuterium evolution by itself is straight forward since D is only produced 
in BBN and destroyed by stellar processing (Reeves et al. 1973; Epstein et al. 1976). Thus 
the present HST determined ISM D/H values give a robust upper bound on rj, but alone 
do not specify the degree of depletion between BBN and the present day. While progress 
has been made recently in determining this primordial value by observing D/H in quasar 
absorption systems, a consistent D p has not yet been found. Both a high primordial value of 
D/H around 2 x 1(T 4 (Carswell et al. 1994; Songaila et al. 1994; Rugers & Hogan 1996a,b) 
and low D p of D/H ~ 2.5 x 1(T 5 (Tytler, Fan & Buries 1996; Buries & Tytler 1996) have 
been measured. Both of these measurements have strengths and weaknesses. The high D 
measurement agrees well with best-fit rj for 4 He and 7 Li (Fields et al. 1996) but requires 
destruction factors of ~ 10 for D to reach its present-day observed values. While models 
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exist which can accomplish this (see e.g. Vangioni-Flam & Audouze 1988; Scully & Olive 
1995; Scully et al. 1996a,b), these models prove problematical for 3 He since in traditional low 
mass star models D is converted into 3 He. The low D measurement on the other hand proves 
less of a problem for traditional 3 He evolution but would require large systematic errors on 
the primordial 4 He measurement (Copi, Schramm & Turner 1995a) and a significantly higher 
primordial 7 Li abundance (by a factor of about 3). The possibility that both the high and 
low D/H measurements are correct and indicate a possible baryon inhomogeneity has been 
explored in Copi, Olive & Schramm (1996). 

Since D is converted to 3 He in the pre-main sequence phase of stellar evolution the 
chemical evolution of these two isotopes is closely linked. In what we refer to as traditional 
models, 3 He may be produced or destroyed, though not totally (no more than what is 
expected in massive stars) in the main sequence phase. As the primordial D/H ratio is 
increased, these models will produce more 3 He, whose evolution must be fit to match the 
solar and present day abundances. We discuss these abundances in §3 below. 

Vangioni-Flam et al. (1994) have explored the evolution of D and 3 He using simple closed- 
box chemical evolution models adopting a D p of D/H = 7.5 x 10~ 5 which is consistent with 
the primordial 7 Li inferred by observations and within 2a plus estimated systematics of 4 He. 
They found that 3 He is overproduced compared with the observed solar and present-day 
values of 3 He unless it is assumed that 3 He is destroyed significantly in low mass stars (i.e., 
at levels comparable to the destruction in more massive stars). There exists, however, some 
observational evidence that at least some low mass stars are net producers of 3 He. Rood, 
Bania, & Wilson (1992) and Rood et al. (1995) find planetary nebulae with abundances of 
3 He as high as ~ 10~ 3 . When 3 He production in low mass stars as calculated by Iben & 
Truran (1978), Vassiladis & Wood (1993), or Weiss, Wagenhuber, & Denissenkov (1995), 
there is always a problem with the overproduction of 3 He and in particular the solar value 
of 3 He/H (Olive et al. 1995; Galli et al. 1995; and Dearborn, Steigman, & Tosi 1996). 

As an attempt to overcome the conflict between these observations and the results of 
Vangioni-Flam et al. (1994), Scully et al. (1996a) proposed that the solution to the prob- 
lem may be one of chemical evolution. They considered models which included a higher 
production of massive stars early in galactic history. This was accomplished by choosing an 
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initial mass function (IMF) which was skewed more towards massive star production at early 
times in galactic evolution but resembled a more normal IMF at later times. Since larger 
mass stars are net destroyers of 3 He this could reduce the predicted abundance at the solar 
epoch. While significantly lower solar 3 He abundances resulted from these models, they still 
exceeded the observed solar value by a factor of ~ 2. They also considered the possibility 
that the solar system 3 He may have been depleted from the galactic abundance at that time 
and computed the degree to which explosions of supernovae in the solar neighborhood may 
have affected the observed elemental abundances. They found that at best only 10% of the 
initial 3 He could be destroyed by this process. They concluded that chemical evolution could 
not fully account for the 3 He problem. 

In an attempt to find models which can reproduce the observed solar and present-day 
D abundances assuming primordial deuterium values as high as those inferred from the 
quasar absorption systems, Scully et al. (1996b) considered a class of models which included 
an early epoch of star formation skewed toward massive star production coupled with a 
supernova wind-driven outflow mechanism. Outflow aids in the D destruction and helps in 
avoiding overproducing metals associated with the massive star production (Vangioni-Flam 
& Casse 1995). These models successfully reproduce a number of observational constraints 
including the age-metallicity relationship and the distribution of low metallicity G-dwarfs. 
These models, however, were found to overproduce 3 He by a factor of ~ 10 at the solar 
epoch assuming standard stellar yields such as those of Iben & Truran (1978). The recent 
observations of Gloeckler & Geiss indicate that the sum (D + 3 He) /H has remained relatively 
constant since the solar epoch. This would indicate that overall low mass stars are neither 
destroyers nor producers of 3 He (Turner et al. 1996). Scully et al. (1996b) found that if 
stars in the narrow mass range of .9 M Q to ~ 1 M Q produce 3 He according to standard 
stellar evolutionary theory while those more massive destroy it, the predicted 3 He abundance 
matches the solar and present-day observations. Furthermore, this would explain the high 
3 He abundances measured in planetary nebulae if their progenitors lie in this mass range. 

Recent work in stellar nucleosynthesis has indicated that low mass stars may indeed be 
net destroyers of 3 He if one includes the effects of rotational mixing in low mass stars on the 
red giant branch (Charbonnel 1994, 1995, Hogan 1995, Wasserburg, Boothroyd, & Sackman, 
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I.- J. 1995, Weiss et al. 1995, Boothroyd & Sackman 1995, Boothroyd & Malaney 1995). In 
contrast to the result of Scully et al. 1996b, the extra mixing does not take place for stars 
which do not undergo a helium flash (i.e. stars > 1.7-2 M Q ). Thus stars less than 1.7 M 
are responsible for the 3 He destruction. This would imply that the progenitor stars of the 
high 3 He observed by Rood et al. must exceed this mass to avoid a conflict. Unfortunately, 
the very high 3 He/H abundances seen in planetary nebulae are more representative of 1 M & 
stars, as these are the only stars able to produce abundances as high as 3 He/H ~ 1CT 3 . (The 
predicted yields of 3 He/H in low mass stars falls as M -2 in standard stellar models (Iben & 
Truran 1978), see also Galli et al. (1996).) 

In this paper, we consider in more detail ways in which low mass stars affect the problem 
of the overproduction of 3 He. We will incorporate yields from stellar models which include 
the effects of rotational mixing (Boothroyd & Malaney 1995) into a wide range of chemical 
evolution models to test their ability to solve the 3 He problem. Note that in these models, 
the destruction of 3 He is not complete. To fit the C and O ratios, only about 80% destruction 
is required after the initial production phase. Thus these "destruction" models actually just 
keep the 3 He approximately constant which fits the observations of Gloeckler & Geiss (1996). 
As an attempt to side-step the apparent conflict between these new yields and the planetary 
nebulae abundances of 3 He, we shall consider the possibility in which some percentage of 
stars (on the order of 10%) follow a post main-sequence evolutionary path which does not 
include rotational mixing. The planetary nebulae observed by Rood et al. would then be 
included in this class of stars. We shall test the ability of this scenario to reproduce the solar 
and present-day 3 He observations. 

Another possibility which we will examine is that the gas expelled from low mass stars 
does not instantaneously mix in with the ISM. The low velocity of the planetary nebulae 
may translate into mixing times of several billions of years. Thus only a fraction of the 
3 He predicted by standard stellar evolution models for low mass stars would be present in 
the gas from which the solar system formed. As we will see, such an assumption is not 
compatible with significant amounts of D destruction. We begin by describing our basic 
chemical evolution models. We briefly outline the observational data of D and 3 He in 
section 3. In section 4, we discuss in detail the results of models which include a time delay 
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in allowing material from planetary nebulae to mix in with the ISM. In section 5 we discuss 
the results of models which include the yields of Boothroyd & Malaney (1995). 

2 Chemical Evolution 

In this section, we provide a basic outline for the chemical evolution models we will be using 
in this work. We want to consider models capable of destroying deuterium by a large factor 
but do not overproduce metals which generally result from the increased stellar processing 
necessary. Scully et al. (1996b) have developed a class of models which include supernova 
wind-driven outflow that satisfy these criteria. In addition, these models satisfy a number of 
observational constraints such as solving the well known G-dwarf problem and reproducing 
the age-metallicity relationship. We shall adapt these models for this work. A summary of 
the main features of these models is given below. Variations about such models such as in 
the stochastic models of Copi, Schramm, & Turner (1995b) may explain the variations in 
3 He in Galactic H n regions but do not change the core results described here. 

To follow the evolution of the mass contained in gas, we use the classical chemical evo- 
lution equations (see e.g. Tinsley 1980). The evolution of the gas mass density in the disk 
including outflow is given by, 

d -^ = -m+e{t)-o{t). (1) 

In this equation, ip(t) is the rate at which gas is being used up by star formation, e(t) is the 
rate at which gas is returned to the ISM by stellar deaths either in supernova events or in 
planetary nebulae given by, 

rm up 

e (^) = / (m — m^^(m,t)^(t — r(m) — T D )dm. (2) 

Here, m up is the upper mass limit of stars that form and m R is the remnant mass. We 
adopt the remnant mass from Iben & Tutukov (1984). r(m) is the main sequence lifetime 
of stars adopted from Scalo (1986). We have also incorporated an additional time delay, T D , 
to account for the delay of material from planetary nebulae in mixing in with ISM. T D is a 
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free parameter in our models for stars with mass less than 8 M and set equal to zero for 
stars more massive than this which we assume will supernova and mix in with the ISM on 
a very short time scale. 

o(t) in equation (H) is the rate at which gas leaves the disk. In this work, we shall adopt 
the outflow mechanism detailed in Scully et al. (1996b). In theses models, some fraction of 
the energy of a supernova event is assumed to go into heating ISM gas to the escape velocity 
of the disk, leaving the system. The rate at which mass is lost from the system can be 
determined from, 

-M w v 2 esc = eE SN N SN , (3) 

where v esc is the escape velocity which is taken to be roughly twice the rotational velocity 
i.e. ~ 500 kms -1 . We assume that the dark matter dominates the gravitational potential of 
the Galaxy and therefore the escape velocity will be taken to be a constant. This will not be 
the case in merger models (see eg. Mathews & Schramm (1993)) where the escape velocity 
would be lower initially allowing for possibly higher outflow rates. E$n is the energy per 
supernova event which is assumed to be 10 51 erg and N S n is the supernova rate. 

All of the supernova energy would be available for heating ISM energy to escape velocity if 
radiative cooling of the expanding shell wasn't important before remnants collide. A critical 
supernova rate can be determined which the actual supernova rate must exceed if cooling 
is unimportant before SNRs collide and merge. David et al. has determined this critical 
supernova rate to be, 

N crit = 0.83kpc- 3 yr \ -) 182 , (4) 

where n is the number density of the ISM gas. Scully et al. (1996b) found that in their 
models Nsn never exceeds this critical value so cooling is always important. Thus in order 
to determine the fraction of energy per event which goes into heating the gas (i.e. e) the 
residual thermal energy of an SNR after it has collided and merged with other remnants 
must be determined. Larson (1974) has estimated the ratio of the residual energy to initial 
supernova energy in terms of the supernova and critical supernova rates to be, 



which is then used to determine the rate of mass loss in equation ([3]). 

In addition to the gas which leaves the disk due to supernova heating, it is assumed that 
some fraction of the supernovae ejecta does not cool radiatively and is flushed directly out 
of the galaxy. Vader (1986) demonstrated that simple supernova-driven wind models with 
a homogeneous ISM cannot reproduce the observed chemical properties of dwarf elliptical 
galaxies and proposed this additional metal enhancement to galactic wind models. The 
fraction of the metals produced in the supernova progenitors which does not cool radiatively 
and is blown out of the galaxy is denoted by v and will be adjusted to match the observed 
metallicity in the solar neighborhood. 

We now turn to the evolution of heavy elements in models including outflow of the type 
described above. We assume that all stars of a mass greater than 8 M & will supernova. It 
will be convenient to rewrite the rate of mass ejected by stars, e(t), to be the sum of the rate 
ejected by those that supernova, e s , and all other stars, e{t) — e s . Equation ([[]) may then be 
rewritten as, 

dMr 

' -i/>{t) + e(t)-ue,(t)-M w , (6) 



dt 

where o(t) has now been replaced by ve s {t) + Mw to reflect the contribution from the wind- 
driven outflow and metal enrichment. Equation (||) may then be extended, 

d ^ G) = + ex(t) - ve sX (t) - M W X, (7) 

where exit) and e s xit) represent the amount of metals ejected by stars and by type II 
supernova respectively. This equation can be further simplified to read, 

dX (ye a (t) - e(t))X - ve aX (t) + e x (t) fo . 



dt M, 



G 



As in Scully et al. (1996b) we will consider three different initial values for primordial 
D/H. In cases in which a high primordial D value is assumed, D/H = 2 x 10~ 4 we will 
consider models similar to model II in Scully et al. (1996b). The evolution of these models 
is divided into two phases. The first phase consists of a steeply declining exponential SFR of 
only massive stars (> 2 M & ) which lasts for only the first .5 - 1.0 Gyr. This is followed by a 
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phase with an exponential SFR with a more normal IMF which continues to the present-day. 
For the first phase, a SFR of the form ij)(t) ~ e - */ T is chosen which lasts for r Gyr. A power 
law IMF is chosen of the form tp(m) ~ m~ 2 " 7 for the range 2 - 100 M . The second phase 
with a SFR ■?/>(£) ~ e - */ 2 ' 5 using the same power law IMF but is in the range .4 - 100 M Q . 
For cases in which an intermediate D p of D/H = 7.5 x 10 -5 is chosen, we shall adopt a 
model similar to model lb of Scully et al. (1996). A SFR proportional to the gas mass is 
chosen, specifically, ip = .28Mg, with a power law IMF, <p oc m" 2 ' 7 in the mass range of .4 - 
100 M . Finally, we will also consider the possibility that the primordial D/H value is low, 
D/H = 2.5 x 10~ 5 . Here we will use the model Ic from Scully et al. (1996b). In this case, a 
constant SFR is assumed ip = 0.07 with the same IMF as in the intermediate D p case. The 
IMF has been normalized in each case such that, 

rrriup 

/ m(p(m)dm = 1, (9) 

where mi ow and m up correspond to the limits of the mass range for each phase. 

We determine the abundance of 16 in our models adopting the yields of Woosley & 
Weaver (1995). As previously mentioned, we have allowed the fraction of supernova which 
participate in the galactic wind, denoted by u, to be a free parameter in all of our models, and 
is adjusted to reproduce the solar abundance of 16 0. In Figure 1, we show the evolution of 
the oxygen abundance corresponding to the models chosen for the three values of (D/H) p and 
are labeled by this value in units of 10~ 5 . Also shown is the evolution of oxygen for the case 
of high primordial D/H in the absence of the winds {v — e — 0) for comparison. Since the 
metal enhanced wind we are employing leads to an enrichment of the intergalactic medium, 
there is a limit to the value of v which can be translated into a limit on the maximum value 
for (D/H) p in these types of models. We will return to this question in §5. All of the models 
we present unless otherwise noted result in a D evolution consistent with the observations 
which we present in the next section. 
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3 Observational Constraints on D and 3 He Evolution 



Since our primary constraint on models of galactic chemical evolution is the abundance of 
D/H (and to some extent 3 He/H), we briefly describe the adopted abundances at the present, 
solar, and in the case of D/H, primordial times. The present day ISM D abundance has been 
determined by Linsky et al. (1993, 1995) using the HST to be 

(D/H) = 1.60 ± .09+^0 x 10 ~ 5 - ( 10 ) 

We shall adopt this value for D/H today. The present day 3 He abundance has been deter- 
mined in a number of galactic Hn regions (Balser et al. 1994). A large range of values exist 
from 1.1 — 4.5 x 10~ 5 . This range may be indicative of a bias and/or pollution (Olive et 
al. 1995) or stochastic evolution (Copi et al. 1995b). A recent measurement of ISM 3 He gives 
(Gloeckler & Geiss 1996) 

( 3 He/H) c = 2.1±;1 x 1(T 5 (11) 

We will use this latter value when comparing to the results of chemical evolution models. 
The presolar D and 3 He abundances were recently discussed in Geiss (1992) and in Scully et 
al. (1996). Our adopted presolar values of D, 3 He, and D + 3 He are: 

[(D + 3 He)/H] & = (4.1 ± 0.6 ± 1.4) x 10" 5 , (12) 

( 3 He/H) & = (1.5 ± 0.2 ± 0.3) x 10" 5 , (13) 

(D/H) Q = (2.6 ± 0.6 ± 1.4) x 10~ 5 . (14) 

We note that recent measurements of surface abundances on Jupiter show a somewhat higher 
value for D/H, D/H = 5±2 x 10~ 5 (Niemann et al. 1996). This value is marginally consistent 
with the inferred meteoritic values. 

Finally, there have been several recent reported measurements of D/H is high redshift 
quasar absorption systems. Such measurements are in principle capable of determining the 
primordial value for D/H and hence r] because of the strong yet monotonic dependence 
of D/H on rj. However, at present, detections of D/H using quasar absorption systems 
indicate both a high and low value of D/H. As such, we caution that these values may not 
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turn out to represent the true primordial value. Nevertheless, we will explore in this work, 
the consequences of choosing either a high or low (D/H) p . Thus we will consider in turn 
(D/H) p = 2.0 x 10~ 4 (Carswell et al. 1994; Songaila et al. 1994; Rugers & Hogan 1996a,b) 
and (D/H) p = 2.5 x 10~ 5 (Tytler, Fan & Buries 1996; Buries & Tytler 1996). We also note 
that the Jovian measurements of D/H are only marginally consistent with the low (D/H) p 
values. 

The evolution of D in the ISM can be determined by extending equation (0). Since D is 
only destroyed in stars, equation fl7|) becomes 

^ = -mo - M W D (15) 

By substitution and some further simplification, we find 

dD = (ue a {t) - e{t))D 

dt Mq 1 ' 

4 Mixing from Planetary Nebulae 

In this section, we test the possibility that 3 He is in fact produced in low mass stars, yet the 
timescale for the mixing of the 3 He enriched gas into the ISM is very slow, on the order of 
several Gyr. In the limiting case, we would withhold all of the gas ejected by stars of under 
~ 8 M & from mixing back in with the ISM. We are motivated to consider delaying this gas 
from mixing in with the ISM since there is in fact some evidence for a dispersion in 3 He in 
galactic Hn regions (Balser et al. 1994). In addition, the long timescales and relatively low 
matter velocity associated with planetary nebulae could contribute to a delay in the mixing 
of 3 He into the ISM. 

To model this effect, we first choose a primordial value for D/H. We then run a chemical 
evolution model such that we delay the gas from stars under 8 M Q from returning to the 
ISM. Figure 2 illustrates the results for high D/H with time delays of 1, 5, and oo Gyr. 
While the 3 He evolution is now improved in each case, none is able to reproduce a viable 
D and 3 He evolution, using traditional assumptions about 3 He evolution in low mass stars. 
Delay times of < 7.0 Gyrs can destroy enough D but still overproduce 3 He. Longer delay 
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times more closely reproduce the 3 He observations but do not destroy D by more than a 
factor of ~ 1.5. 

Figure 3 shows similar results for model I with the lower (D/H) p = 7.5 x 10~ 5 . While a 
delay time of oo can reproduce the 3 He observations, again D can not be destroyed by more 
than a factor of ~ 1.5. We estimate that the longest time delay which can still reproduce 
both the solar and present-day D observations can only reduce the 3 He/H by a factor of 
20% from the same model run with no time delay and traditional 3 He evolution in low mass 
stars. We can conclude that this effect can clearly not be entirely responsible for the observed 
solar value of 3 He. This combined with our previous results (Scully et al. 1996a and Scully 
et al. 1996b) all point to the stellar yields as being primarily responsible for the observed 
3 He/H abundances. 

The exception to this conclusion is the case of low primordial D/H. Of course in this case, 
D need not be destroyed by much more than a factor of 1.5. In Scully et al. (1996b), this 
case did not greatly overproduce 3 He, and a time delay of even 1 Gyr would further improve 
the comparison with the data. It remains a difficulty in these models however, to produce 
sufficient amount of heavy elements such as oxygen. 

5 Effects of New 3 He Yields 

A number of recent papers suggest that stars less massive than 2 M may be net destroyers 
of 3 He when the effects of rotational mixing in the red giant phase of stellar evolution is 
included (see e.g. Charbonnel 1994, 1995, Hogan 1995, Wasserburg, Boothroyd, & Sackman, 
I.- J. 1995, Weiss et al. 1995, Boothroyd & Sackman 1995, Boothroyd k Malaney 1995). 
Boothroyd & Malaney (1995) give detailed results for 3 He destruction factors which we have 
implemented into our chemical evolution models. Main sequence mass loss is not expected 
to affect these yields due to the small mass loss rates for low mass stars. We remind the 
reader that these 3 He destruction mechanisms take place in the post-main sequence phase 
of stellar evolution. Traditional main sequence evolution is unchanged. 

Boothroyd & Malaney have already tested their results using the models of Vangioni- 
Flam et al. 1994 and have adopted a model in which ip = .25a. They have suggested that 
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values of (D + 3 He) p higher than 1.2 x 10~ 4 are unable to reproduce a viable D + 3 He 
evolution using their 3 He yields. This corresponds to an upper limit to (D/H) p < 1.0 x 10~ 4 . 
The model they have chosen, however, is designed to give the correct gas fraction and D 
destruction for a (D/H) p = 7.5 x 10~ 5 . In the following, we will examine the evolution of D 
and 3 He in models which are more suited to destroy D for the chosen D p as well as satisfy 
other observational constraints including the present-day gas fraction using the Boothroyd 
& Malaney yields for 3 He. 

In Figure 4, we show the results of including the Boothroyd & Malaney yields into model 
II of Scully et al. (1996b) with high primordial D/H. As one can see the model has been chosen 
to match the solar and ISM D/H abundances starting at (D/H) p = 2 x 10~ 4 . As we indicated 
earlier, this model incorporates supernovae driven winds to regulate the oxygen abundance. 
The evolution of 3 He/H using the Boothroyd & Malaney yields is shown by the solid curve in 
Figure 4. This model is able to reproduce the present day 3 He/H observation and is within 
the errors including systematics for the solar observation. The effects of rotational mixing 
on the red giant branch on the 3 He/H abundance is dramatic as is seen by comparing the 
solid curve to the dashed one, which utilizes the standard stellar yields of Iben & Truran 
(1978). 

We have also tested the new yields for the case of intermediate (D/H) p = 7.5 x 10~ 5 . 
Figure 5 shows the resulting D/H and 3 He/H evolution. In this case, the solar 3 He/H 
observation is reproduced but the model gives a somewhat low abundance of 3 He/H with 
respect to the present day observation. Once again, we see the significant effect of the yields 
by comparing the solid and dashed curves for 3 He/H. For completeness, we show the results 
for the case of low (D/H) p = 2.5 x 10~ 5 in Figure 6. In this case, because low primordial D/H 
corresponds to both low rj and low primordial 3 He/H, the Boothroyd and Malaney yields 
predict too little 3 He/H. Indeed, the standard yields as shown by the dashed curve show an 
evolution closer to the data. Curiously, it seems that the reduced yields in fact work better 
for the higher Dp. We can conclude that the new yields based on 3 He destruction on the 
Red Giant Branch of low mass stars can solve the 3 He problem even (and especially) for the 
higher D p , when combined with the Galactic evolution models of the type discussed. 

Although the reduced yields of Boothroyd & Malaney (1995) are capable of explaining 
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the rather flat evolution over time of 3 He/H, they can not account for the high 3 He/H 
observed in planetary nebulae (Rood et al. 1992,1994). As Boothroyd & Malaney suggest, 
it may be that not all stars undergo the extra mixing subsequent to first dredge-up and 
that the main-sequence produced 3 He remains intact. In Scully et al. (1996b), it was shown 
that the 3 He evolution and planetary nebulae data could be fit if only stars between 0.9 and 
~ 1M & produced 3 He in significant quantities. For the initial mass functions considered, this 
represented some 10 - 15 % of all stars becoming planetary nebulae. In Figure 7, we show 
the result for the high D/H case when 90% of all stars in the 1-8 M Q mass range undergo 
the post main-sequence destruction while 10% do not. It is this 10% which accounts for the 
observation of high 3 He in planetary nebulae. Recall that all stars in this mass range produce 
3 He while on the main sequence. In most cases, (chosen here to be 90%) post-main sequence 
processing leads to the destruction of 3 He. In the remaining cases, 3 He destruction might be 
stopped is some binary systems or in systems in which the red giant envelope configuration 
is disturbed. The dashed line in Figure 7 shows the case when 100% of the stars destroy 
3 He as given by the reduced yields. As one can see, the reduced yields seem to fit better, 
slightly overproducing 3 He at the solar epoch and slightly underproducing 3 He today. With 
the 10% mix of stars which do not destroy 3 He, 3 He is overproduced at the solar epoch by 
a factor of about 2, though the present day abundance is acceptable. 

In Figure 8, we show the cases for (D/H) p = 7.5 x 10~ 5 with a 90%/10% mix as discussed 
above and in Figure 9, for the case of low D/H, the mix is 40%/60%. Again, the dashed 
curves show the results using only the reduced yields. As one can see from these figures, the 
addition of some stars which do not undergo the 3 He destruction processes on the red giant 
branch are quite consistent with the solar and ISM observations of 3 He. 

Finally, as we noted earlier, a considerable amount of enriched material is expelled from 
the galaxy in these types of models, particularly those with high primordial D/H. Such a 
metal enhanced wind would contribute to the enrichment of the extra-galactic gas with heavy 
elements and indeed there is evidence that such an enrichment has occurred. For example, 
in the clusters observed by Mushotzky et al. (1996) and Loewenstein & Mushotzky (1996) 
the mean oxygen abundance was found to be roughly half solar. For our model with high 
(D/H) p , the enrichment was found to produce a metallicity of about Zi GM ~ 0.2Z Q (Scully 
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et al. 1996b) assuming a total intergalactic gas mass equal to ten times the Galactic mass. 
We have explored increasing (D/H) p in our models to determine the maximum allowed value 
without overproducing metals in the intergalactic medium due to our outflow mechanism. 
For (D/H) p = 3 x 10~ 4 , Zigm ~ 0.3Z & . In our models, we can not increase the deuterium 
abundance beyond (D/H) p = 5.4 x 10~ 4 as higher primordial abundances would require 
such an increased SFR that no gas would remain in the disk today. This value of (D/H) p 
corresponds to a Ziqm ~ 0.6Z Q . 

In completing this paper, we became aware of the very recent work of Galli et al. (1996). 
In that work, Galli et al. (1996) determine the masses of the planetary nebulae observed by 
Rood et al. (1992, 1995). Indeed, they confirm that the progenitor masses of the nebulae are 
<, 2M , as was suspected. They also consider a chemical evolution model with (D/H) p ~ 
3.5 x 10~ 5 with standard stellar yields as well as those provided by Boothroyd which account 
for 3 He destruction on the red giant branch. To match the solar and present day abundances 
of 3 He, they too require a mixture of these yields and their results are in qualitative agreement 
with those presented here. 

6 Conclusion 

It is clear that the evolution of 3 He is complicated by our uncertainties in both the galactic 
and chemical evolution of this isotope. It appears that the abundance of 3 He is in fact little 
changed over the history of the Galaxy. This is particularly difficult to understand if the 
primordial abundance of deuterium is as high as observed in some recent measurements of 
D/H in quasar absorption systems because chemical evolution models using standard stellar 
yields for 3 He in low mass stars are known to lead to a gross overproduction of 3 He at both 
the solar and present day epochs. Unless the total deuterium astration factor is less than 
about 1.5, we have found that imposing a delay for the mixing of 3 He rich material does 
not substantially improve the evolution of 3 He. This further reinforces our conclusion that 
chemical evolution alone can not solve the 3 He problem. 

In this paper, we have implemented new 3 He yields for low mass stars which account for 
possible extra-mixing mechanisms on the red giant branch and lead to a strong depletion 
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of 3 He (Boothroyd & Malaney 1996). These yields do in fact provide for a flat evolution 
for 3 He over the history of the Galaxy. However, they can not account for the high 3 He 
abundances observed in planetary nebulae. As such, we argue that although these extra- 
mixing mechanisms may be operative for most low mass stars in their post-main sequence 
evolution, in some fraction of stars the main sequence produced 3 He must remain in tact. 
In some cases, such a mix, fits the data (solar and ISM) quite well. We also determine an 
upper bound to (D/H) p < 5x 1CT 4 based on the maximum allowed amount of extra-galactic 
heavy element enrichment. This limit assumes that the extra-galactic medium has less than 
one half solar metallicity. This limit is significantly higher than the highest values of D/H 
observed in quasar absorption systems. In conclusion, the present day and solar abundances 
of D and 3 He can be made consistent with either high or low primordial D/H (low or high 
if). Thus the solution to the the " 3 He problem" is clearly not cosmological but rather of 
Galactic or stellar origin. 
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Figure Captions 



Figure 1: The evolution of 16 O/ 16 O for the chemical evolution models described in the 
text corresponding to (D/H) p = 2 x 10~ 4 with outflow (solid line) and without outflow 
(dotted line), (D/H) p = 7.5 x 1(T 5 (dashed line) and (D/H) p = 2.5 x i(T 5 (dot-dashed 
line). 

Figure 2: The evolution of D/H and 3 He/H assuming a high primordial D/H of 2 xlO~ 4 . 
Solid curves show the results of standard model 3 He yields and the chemical evolution 
model (II) of Scully et al. (1996b). Shown in this figure is the resulting evolution when 
a delay is applied to the mixing time for the gas ejected by planetary nebulae. Delays of 
1, 5, oo are shown by the dotted, dashed, and dot-dashed curves respectively. 

Figure 3: As in Figure 2 for (D/H) p = 7.5 x 1CT 5 . 

Figure 4: The evolution of D/H and 3 He/H assuming a high primordial D/H abundance 
of 2 xlO~ 4 . Solid curves show the results of the reduced 3 He yields from Boothroyd 
& Malaney (1995) and the chemical evolution model (II) of Scully et al. (1996b). The 
dashed curve shows the result using standard 3 He yields. 

Figure 5: As in Figure 4 for (D/H) p = 7.5 x 1CT 5 . 

Figure 6: As in Figure 4 for (D/H)„ = 2.5 x 1CT 5 . 

Figure 7: The evolution of D/H and 3 He/H assuming a high primordial D/H of 2 xlO" 4 . 
Solid curves show the results of assuming that 90% of the stars in the 1-8 M range 
destroy 3 He as in the reduced 3 He yields of Boothroyd & Malaney (1995), the remaining 
10% of stars in this mass range produce 3 He using standard stellar yields. The dashed 
curve shows the result from (the solid curve of) Figure 2. 

Figure 8: As in Figure 7 for (D/H) p = 7.5 x 10' 5 . 

Figure 9: As in Figure 7 for (D/H) p = 2.5 x 10~ 5 . Here, the mixture is 40% reduced yields, 
60% standard yields. 
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